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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
ADVANCE RESTRICTED REPORT 

AN INVESTIGATION OF A THERMAL ICE-PREVENTION SYSTEM 
FOR A C-46 CARGO AIRPLANE 
I - ANALYSIS OF THE THERMAL DESIGN 
FOR WINGS, EMPENNAGE, AND WINDSHIELD. 

By Carr B. Neel, Jr. 

SUMMARY 


This report, the first of a series on an investigation 
to develop a thermal i c e- prevent i on system for the C-46 cargo 
airplane and to provide fundamental information concerning 
the performance and operation of the system, presents an 
analysis of the design cf the thermal anti-icing equipment 
for the wings, the empennage, and the windshield. The pur- 
pose of the analysis was to establish the thermal require- 
ments for the C-46 ice-prevention equipment on the basis of 
a certain temperature rise for all points on the heated sur- 
faces of the wings and the empennage, and of a specified heat 
flow through the outer panel of the windshield. The rate of 
■"heat transfer from the wing and empennage surfaces was calcu- 
lated by utilizing the airfoil pressure distributions for 
these surfaces measured in flight. 

The design analysis showed that the thermal requirements 
for the anti-icing systems of the wings, the empennage, and 
the windshield were within the limits of practicability, al- 
though complete protection of the windshield might not be re- 
alized in the specific application when utilizing a secondary 
heat exchanger . 


INTRODUCTION 


For several years the NACA has been engaged in a re- 
search program to investigate the feasibility of utilizing 
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the waste heat of airplane engine exhaust gases to heat those 
parts of an airplane that require protection from the forma- 
tion of ice in order to provide safe and efficient operation 
of the airplane in icing conditions. In this thermal anti- 
icing method, the heat of the exhaust gases is transmitted to 
air in suitable heat exchangers located adjacent to the 
engines t and the heated air is then circulated through the 
airplane so that it passes along the inner surfaces of af- 
fected parts. The method has been demonstrated to be sound 
in successful applications to a Lockheed 12A airplane, a 
Consolidated B-24D airplane, and a Boeing B-17F airplane 
(references 1, 2, and 3, respectively). These applications 
were made at the Committee's Ames Aeronautical Laboratory, 
Moffett 'Field, Calif. 

As a continuation of this general research program on 
aircraft ice prevention, the laboratory has undertaken a com- 
prehensive investigation to develop a thermal ice-prevention 
system for a Curt i ss-Wri ght C-46 cargo airplane. The inves- 
tigation consists of a design analysis of the heat anti-icing 
system, design and construction of the exhaust-gas-t o-air 
heat exchangers, fabrication and instrumentation of the com- 
plete system, and performance tests of the C-46 airplane 
equipped with the thermal ice-prevention system under dry air 
and natural icing conditions. 

This report is the first of a series describing the in- 
vest i gat i on , and presents the design analysis of the thermal 
ice-prevention system for the wings, the empennage, and the 
windshield. Emphasis was placed on this phase of the inves- 
tigation because it was believed necessary to establish a 
satisfactory theoretical basis for the design of the system 
before extensive application of the system could be accom- 
plished. The investigation has been undertaken at the re- 
quest and with the cooperation of the Air Technical Service 
Command of the U. S. Army Air Forces and with the assistance 
of the Curt i ss-Wright Corporation. 


TENTATIVE DESIGN OF THE THERMAL ICE-PREVENTION SYSTEM 


A standard C-46 cargo airplane, shown in figure 1, was 
provided for the investigation. The tentative design of the 
thermal ice-prevention system, on which the subsequent anal- 
ysis was made, was based on previous experience (references 
1, 2, and 3). The important elements of the design are shown 
in figures 2, 3, 4, 5, and 6. 'The existing structure of the 
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airplane and other practical considerations limited various 
factors of the design. 

/ 

In this design the required supply of heated air for 
the ice-prevention system was provided by four exhaust-gas- 
to-air heat exchangers, two in each engine nacelle. The 
heat exchangers were located on the inboard and outboard 
sides of the nacelles so that the two exhaust-gas pipes from 
the collector rings on the engines led to the exhaust -gas 
entrances of the heat exchangers. On the basis of the design 
the outside air was passed through the exchangers and thence 
circulated through ducts to the wings, the empennage, and the 
windshield. The exchangers on the outboard sides of the na- 
celles supplied heated air to the wing outer panels, and the 
two inboard exchangers supplied heated air to the wing inooard 
panels, the empennage, and the windshield. The disposition 
of the heated air as considered in the analysis is shown in 
figure .2. 

The proposed construction of the outboard wing leading- 
edge regions, which in this report include the airfoil sec- 
tions forward of the 10-percent chord points, is illustrated 
in figure 3. Available dies made the corrugation-type passage 
shown in figure 3 most expedient for the internal circulation 
of heated air over the surfaces to bp protected. The upper 
and lower corrugations were divided by a continuous spanwise 
slot at the 0-percent reference chord point in order to allow 
the entrance of heated air, whi ch flowed chordwise between 
the corrugations and the outer wing skin. The location of 
essential wing structure at 10-percent chord required that 
the corrugations be terminated at that point. The baffle 
plate, which was located between the 5- and 6-percent chord 
points (because this location was most advantageous from the 
standpoint of fabrication) formed, with the leading-edge 
inner surface, the spanwise duct to distribute heated air to 
the corrugations. , 

Structural considerations required t hat most of the nose 
ribs be retained. However, the installation of a corrugated 
inner skin is known to add a certain amount of rigidity to 
wing and empennage structures, and accordingly the nose ribs 
outboard of station 292 (fig. 2) were omitted. A liner was 
installed along the inner edges of the nose ribs in order to 
give a smoother contour to the wing leading-edge duct than 
would be possible if the nose ribs were exposed to offer re- 
sistance to the flow of heated air. The liner was terminated 
at station 292 because the design called for omission of the 
nose ribs outboard of that point. Holes were cut in the nose- 
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rib liner as shown in figure 3 to allow the passage of heated 
air from the leading-edge duct to the corrugation entrances. 
Thus the heated air passed along the leading-edge duct (re- 
gion 1 of fig. 3), through the corrugation gaps (region 2) 
into the space between the baffle and front spar (region 3), 
through the holes in the spars, and from the afterbody of the 
wing (region 4) out through exit holes in the aileron and 
flap-hinge regions (region 5). 

The distribution of heated air for the empennage and de- 
tails of the leading-edge region construction of these sur- 
faces, as considered in the analysis, are shown in figures 4 
and 5, respectively. The proposed nor i zontal- stabili zer and 
vertical-fin thermal systems are similar to that of the wing 
outer panels. However, the nose ribs were eliminated entirely 
and the leading-edge ducts were formed by the corrugations and 
the baffle plate.' 

The tentative windshield thermal system design is illus- 
trated in figure 6, which shows the pilot's windshield only. 
This design consists of two glass panels spaced with a con- 
stant gap between them. The existing windshield outer panels 
are made up of two layers of semi-t empered plate glass sepa- 
rated by a layer of vinyl plastic. In the tentative anti- 
icing-system design the heated air entered the plenum chamber 
at the base of the windshield panels, passed upward through 
the gap between the outer and inner glass panes, and exited 
at the top into the cabin, thus heating the pilot's compart- 
ment as well as protecting the windshield from icing. Because 
of this special design feature, however, it was necessary to 
include in the windshield anti-icing system a secondary heat 
exchanger in order to prevent carbon monoxide from entering 
the cabin in the event of a failure of the exhaust -gas-t o-air 
heat exchanger. The secondary heat exchanger transferred 
heat from the primary heating air, supplied by the inboard 
exhaust-ga s-t o-air heat exchangers, to the secondary air from 
the ambient stream. The windshield ice-prevention system in- 
cluded valves by means of which the secondary heated air could 
be discharged directly into the pilot's compartment during 
periods when windshield heating is not necessary. 


TEST PROCEDURE 


The airfoil pressure distributions over the wings and 
the empennage were obtained experimentally in flight in order 
to provide a basis in the analysis for calculating the heat 
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transfer from the heated surfaces of the airplane to the out- 
side air. This method of analysis is believed to be a more 
rational procedure for obtaining the heat requirements of 
these surfaces than the methods used in previous designs of 
aircraft thermal ice-prevention systems. 

The pressure distributions were measured by means of 
pressure belts, as described in reference 4. Data were ob- 
tained for stations 82, 157, 292, and 382 on the outer panel 
of the right wing; stations 72. 127, and 172 on the right 
stabilizer; and stations 122 and 172 on the fin. The airfoil 
sections of the wing consist of an NACA 23017 section of 16,5— 
foot chord at station 0, tapering to an NACA 4410.5 section 
of 7.1-foot chord at station 393, which is 19 inches inboard^ 
of the wing tip joint (station 412). The wing at station 39o 
has three degrees of washout with respect to station 0. The 
horizontal stabilizer and the vertical fin both have NACA 0010 
airfoil sections throughout. 

The pressure belts were wrapped around the leading edges 
of the wing and empennage surfaces and extended back approxi- 
mately to the 20-percent chord point. Typical installations 
are shown in figures 7 and 8, and the water manometer used 
for indicating pressures at the empennage pressure belts is 
shown in figure 9. 

Previous experience with similar thermal ice-prevention 
installations has indicated that a satisfactory assumed de- 
sign flight condition is cruising engine power for maximum 
range at a pressure altitude of 18,000 feet. The first two 
flights were made at a pressure altitude of 18,000 feet with 
engine-power setting to approximate maximum range of the air- 
plane, and the pressure distributions for the two outboard 
stations of the wing and the horizontal stabilizer were ob- 
tained. The indicated airspeed was established at 155 miles 
per hour under level-flight conditions. The required indi- 
cated airspeed having been determined, the remaining pressure 
distributions were secured at 5000 feet pressure altitude and 
approximately 155 miles per hour indicated airspeed. It was 
assumed that at the same indicated airspeed the angle of at- 
tack of the airplane and the airfoil pressure distribution 
are unchanged with variations in altitude for level flight at 
a constant load condition. The gross weight of the airplane 
and crew, in flight, for all the pressure-belt tests was ap- 
proximately 39,000 pounds. During the flights, the water 
manometers indicating the pressures at the pressure belts 
were allowed to reach equilibrium, and the manometers were 
photographed to record the readings. Adter each flight, the 
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pressure belts were moved to the next station and the pro- 
cedure was repeated. 


SYMBOLS 


For the purpose of analysis, the following symbols were 

used: 

A cross-sectional area, square feet 

c length of airfoil chord, feet 

Cp specific heat of air, British thermal units per pound, 
degree Fahrenheit 

D e equivalent diameter of a duct, equal to 4m, feet 

f friction factor for air flow in ducts 

g gravitational acceleration, taken as 33.2 feet per sec- 
ond, sec ond 

G weight rate of air flow per unit of cross-sectional area 
pounds per second, square foot 

h surface heat -transfer coefficient, British thermal units 
per hour, square foot, degree Fahrenheit 

k thermal conductivity of air, British thermal units per 
hour, square foot, degree Fahrenheit per foot 

L distance around airfoil surface, measured chordwise from 
O-percent-chord point, feet 

m hydraulic radius, or ratio of cross-sectional area to 
wetted perimeter in a duct, feet 

N length of air duct or corrugation, feet 

Nu Nusselt number, equal to hD e /k 

p calculated static pressure, pounds per square foot or 
inches of water 

p^ measured pressure differential from free-stream static 
pressure, inches of water 
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p 0 free- stream static pressure, pounds per square foot 
P pressure coefficient, equal to Pf/q 0 
q 0 free-stream dynamic pressure, inches of water 

Q, rate of heat flow, British thermal units per hour 

R gas constant for air, taken as 53.3 feet per degree 

Fahrenheit 

R c "body Reynolds number, equal to cV 0 /y 
Re Reynolds number, equal to G-D e /y> 

s distance around airfoil surface, measured chordwise 
from stagnation point, feet 

S surface area, square feet 

t temperature, degrees Fahrenheit 

T temperature, degrees Fahrenheit absolute 

v specific volume, cubic feet per pound 

V local ambient-air velocity outside the boundary layer, 
feet per second 

V\ airplane indicated airspeed, miles per hour 


V o 

fr ee- st ream 

velocity, feet 

per second 


w 

weight rate 

of 

air flow per corrugation, 

pounds per hour 

w a 

weight rate 

of 

air flow, pounds per 

hour 


X 

distance al 

ong 

chord line, 

measured 

from 

O-percent- chord 


point 






y 

density of 

air , 

, pounds -p er 

cubic fo< 

Dt 


8 

boundary- layer 

thi ckne ss , 

feet 




X constant dependent on shape of boundary-layer velocity 
profile 
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M> absolute viscosity of air, pounds per second, feet 

p kinematic viscosity of air, equal to p/V, square feet 

per second 

For the sake of clarity, other symbols, not presented 
here, will be defined when introduced. 

The values of Cp , |x, and k for air, subsequently 

used in the calculations, were taken from figure 10, which 
contains plots of data obtained from reference 5. 


Use of Subscripts 

The temperature, pressure, and specific volume subscript 
numbers refer to the temperature, the pressure, or the spe- 
cific volume of the region to which the numbers apply (figs. 

3 and 5). Thus, the temperature of the air in a corrugation, 
region 2y will be designated as 1 2 , the static pressure as 
p 2 , and the specific volume as v 2 . 

Two subscript numbers refer to the flow of heat from one 
region to another. Thus, the heat-transfer coefficient for 
the air in a corrugation to the wing surface will be desig- 
nated as h 

2 — 6 

The subscript av refers to average conditions. 


ANALYSIS OF WING AND EMPENNAGE THERMAL SYSTEMS 


Past experience has indicated that a satisfactory design 
specification to insure ice prevention is the maintenance of 
the heated surfaces at appr oximat ely 100° F above ambient-air 
temperature during flight in dry air (no visible moisture) at 
a pressure altitude of 18,000 feet with cruising engine power 
for maximum range. Therefore, for the purpose of this anal- 
ysis, the requirement of a minimum temperature rise of approx- 
imately 100 F for all points on the surface of the wing and the 
empennage between the 0- and 10-percent-chord points was 
established for the design flight conditions of 0° F ambient- 
air temperature, 18,000 feet pressure altitude, and 155 miles 
per hour indicated airspeed. 

Preliminary investigation of the structure in which the 
system was to be installed indicated that essential wing 
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structure precluded freedom of design. An analysis indicated, 
however, that a satisfactory thermal design could he obtained 
without altering this essential primary structure, and use 
could be made of existing die equipment for fabricating the 
inner skin. As a result, certain design factors were consid- 
ered fixed at the start of the analysis, such as the extent 
of the heated region or the length of the corrugations, which 
was fixed by the location of essential structure at the 10- 
percent-chord point, and the corrugation size, which was de- 
termined by available dies. With these factors fixed, the 
analysis resolved itself into the adjustment of such variable 
factors as the heated-air flow rate and the pressure drop to 
satisfy the design specification. The analysis therefore 
consisted of the determination of (l) the boundary-layer 
characteristics from which the external heat-t ran sf er coeffi- 
cients could be calculated, ( 2 ) the adjustment of the heated- 
air flow rate to obtain the required surface temperature 
rise, and (3) the calculation of the resulting air-pressure 
drops to ascertain that equilibrium had been established. 

These phases are discussed hereinafter under appropriate sec- 
tions of the report. 

A trial- and- error process was employed in the calculation 
of the thermal and pressure values. A more elegant analysis 
was precluded by the great number of variables involved and 
the complicated relationships resulting when expressed in 
analytic form. 

The analysis was made for the wing outer panel from sta- 
tion 11 to station 412, for the horizontal stabilizer from 
station 47 to station 190, and for the vertical fin from sta- 
tion 112 to station 192. The analytical study did not in- 
clude consideration of the wing inboard panel or the tips of 
the wing, the stabilizer, or the fin. 


Boundary-Layer Characteristics 

The determination of the external heat-transfer coeffi- 
cients is dependent upon knowledge of the characteristics of 
the boundary layer surrounding the heated surfaces. The 
pressure-belt data recorded during flight tests were used to 
calculate the coefficients of heat transfer from the airfoil 
surface to the boundary layer at each belt location. The 
differences between the measured static pressures at the belt 
orifices and the free-stream static pressure were determined 
and are designated as pf. Table I presents all calculated 

values of (l - P) where P = Pf/q 0 . A plot of (1 --?) as a 
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function of the chord in percent for wing station 82 is 
given in figure 11. Similar plots of airfoil pressure dis- 
tribution were made for each set of pressure-belt data, but 
are not presented in this report. Values of (l - P) were 
taken from these curves for convenient chord stations and 
corresponding values of the velocity ratio-, V/V Q , were 

calculated by the equation 

V/V 0 =/ 1 - P (1) 

which is derived from Bernoulli's equation for fluid flow. 
Figure 12 gives a curve of the velocity ratio thus obtained 
for flow over the wing at station 82. 

The location of the laminar- separat i on point was deter- 
mined as a function of the slopes of the velocity profile. 
The curves of velocity ratio were replaced by two straight 
lines (referred to as double-roof profiles in reference 6) 
the slopes of whi ch closely approximated the slopes of the 
actual velocity profile. The ratios of the slopes of the 
double-roof profiles were calculated, and the ratios of the 
local velocity at separation to the maximum local velocity 
were obtained from figure 13, where b is the slope of in- 
creasing velocity and 8 is the slope of decreasing veloc- 
ity. Figure 13 is a plot of the following equation which 
was obtained from footnote 1 on page 21 of reference 6, 

^ =/l - r (6 - 1) (2) 

U ^ 5 


wher e 

U s local velocity at laminar separation 

U x maximum local velocity 

r ratio of 8 to b 

8 S a function of r 

For the upper surface at wing station 82, 


8 = 0.278 


b 


6.00 
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Then 


r = £ = = 0.04 63 

b 6.0 0 


and from figure 13, 




= 0.906 


From figure 12, 


then 


max 


u i 

— = 1.52 
v o 


U,' = 1.52 V, 


and 


U 


S = 0.906. X 1.52 = 1.378 


Thus, by extending the decreasing velocity curve beyond the 
limits of the graph shown in figure 12, it was found that 
separation occurred at s/c x 100 = 27. 

Transition from laminar to turbulent flow over the air- 
foil surfaces was considered to take place at the laminar- 
separation point. In all cases, the separation point was 
established well aft of the region to be heated. Since 
the boundary layer was found by the calculations to be lam- 
inar over the heated areas, the method presented in reference 
7 was used to determine the boundary-layer thickness. The 
equation for boundary-layer thickness 6 in the laminar-flow 
region, as given in reference 7, is 

r *i/e 1 



where the subscript l denotes the point on the airfoil sur- 
face at which the boundary-layer thickness is being computed. 
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si/c 


8 C 17 


The evaluation of J ^ (^) was made graphically in 

o 0 

each case, using curves of the functions of velocity ratio 
similar to those shown in figure 14. The values of the con- 
stant • 5.3 c /E c were determined from the true airspeed and 

the chord length. 

For wing station 82, 

c = 14.51 feet 


The indicated airspeed deviated slightly from the desired 
value of 155 miles per hour during each flight when the pres- 
sure-belt data were recorded. The indicated airspeed was 151 
miles per hour during the flight when belt data were obtained 
for station 82. Correcting the indicated airspeed to free- 
stream velocity, using the reciprocal square root of density 
ratio at 18,000 feet pressure altitude, gives 

V Q = 1.325 "V ^ = 1.325 x 151 = 200 miles per hour 

= 293,5 feet per second 


At 18,000 feet pressure altitude and 0° F free-air tempera- 
ture, 


fy _ P o 1058 

53.3 X 460 


0.0431 pounds per cubic foot 


From figure 10, 




= 1.104 X 10 pound per second, feet at 0° F 


Then 


,, _M>_ 1.104 X 10 

y - y = 6 ”0431 ~ 0 * 00^256 square foot per second 


R c = 


cV 14.51 X 293.5 


0 . 000256 


— = 16.64 X 10 


and 


5.3c 5 _ 5 .3 X 14 . 51 5 
R c 16.64 X 10 6 


0.0000671 
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Since the value of 8 is indeterminate at the stagna- 
tion point, the following appr oximat i on recommended in refer- 
ence 7 was used: 

8 s I : (4) 

stag 5H C c 


where r is the radius of curvature at the stagnation point 
For wing station 82, 

F = 5.26 inches 


and 


14 .51 


5.26 


stag 


5 X 16.64 x 10' 


14.51 X 12 


= 7.65 X 10 


- e 


6 stag = 2,76 x 10 feet 


The values of r were appr oximated from the airfoil-section 
loft lines. A curve of the calculated "boundary-layer thick- 
ness for wing station 82 is shown in figure 12. 


Calculation of Heat-Transfer Coefficients 


The coefficient of convective heat transfer from a 
heated srirface to the laminar boundary air layer can be ex- 
pressed as a function of the boundary-layer thickness by the 
e quat i on 


h 


6-7 


kX 

6 


(5) 


The value of X was considered to be a constant (0.765) 
from the stagnation region to the minimum pressure point, 
and then vas considered to decrease linearly from 0.765 to 0 
at the laminar- separat i on point. (See reference 7.) 

Since the indicated airspeed varied slightly for each 
flight when pressure-belt data were obtained, the calculated 
values of h c „ were corrected to give corresponding values 

at the design airspeed, 155 miles per hour, by the equation 
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(h 6 _ 7 ) 


cor r 



( 6 ) 


A tabulation of the results of intermediate steps taken to 
calculate the boundary-layer thickness and the external heat- 
transfer coefficient for wing station 82 is presented in 
table II, The curves of calculated values of corrected ex- 
ternal heat-transfer coefficients from the pressure-belt data 
are presented for all wing and empennage stations investi- 
gated in figures 15 to 23. Because the flow was assumed to 
be symmetrical over the vertical fin, the heat-transfer coef- 
ficients were calculated for one side only, as shown in fig- 
ures 22 and 23. The curves sho\* that in all cases the exter- 
nal heat-transfer coefficient increases as the chord and the 
leading-edge radius decrease, with the exception of station 
127 of the horizontal stabilizer, at which point the calcu- 
lated heat-transfer coefficient is highest for the stabilizer 
It is believed that this increased heat-transfer coefficient 
for station 127 of the stabilizer is due to the effect of the 
higher air velocities in this region caused by the propeller 
slipstream . 


Air-Flow and Heat-Flow Calculations 

The weight of air required to produce the desired skin- 
temperature rise was evaluated on the basis of tentatively 
fixed dimensions and coverage of the corrugations and the 
allowable inlet-air temperature. The conditions, which must 
exist for the establishment of a satisfactory flow distribu- 
tion through the system, are that the pressure losses along 
all streamlines must be equal and that the skin-temperature 
rise must provide ice protection. Pressure equilibrium for 
the system was established by calculating the losses in the 
corrugations at the wing-root and wing-tip sections and ad- 
justing the baffle position to provide spanwise duct loss 
equal to the difference between the two chardwise corrugation 
losses. 

For the calculations to determine the required heated- 
air flow to produce a desired temperature rise, the leading- 
edge region was divided into three segments as follows: 0- 

to 3-percent chord, 3- to 7-percent chord, and 7- to 10-per- 
cent chord. Since the corrugation passages were of constant 
cross section, the exact value of the design skin-temperature 
rise (100° F) could not be established at all points on the 
leading edge - that is, some variation was unavoidable. The 
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procedure followed was to assume 100° F rise for the skin 
from 0- to 3-percent chord, calculate the required heated- 
air flow, and then determine the temperature rise aft of 3-. 
percent chord produced “by that air flow. After satisfactory 
air flows were established for the corrugation passages at 
the wind root and tip, the baffle location was determined. 

The establishment of pressure equilibrium, which was accom- 
plished by a trial-and-err or method, required slight changes 
in the heated-air distribution which, in turn, produced 
changes in the skin- temperature distribution. These changes, 
however, were not excessive and the final skin temperatures 
all were appr oximat ely equal to the design requirement. The 
detailed application of the general procedure outlined above 
for the wing outer-panel and empennage leading-edge systems 
is presented in sufficient detail to illustrate the process; 
the remainder of the analysis is presented in tabular form. 

First approximation of required heated-air di str ibut iotu - 
The first calculation consisted in establishing the heated- 
air flow required at two stations, one at the wing root and 
one at the wing tip, to produce a skin-temperature rise of 
100° F from 0- to 3-percent chord. Considering the flow in a 
passage formed by a single corrugation on the upper surface 
at station 82, the heat flow from a 1-inch strip of skin to 
the ambient air must be balanced by an equal flow of heat 
from the air in the passage to the skin, or 

^ 6-7 ® ^ ^ 6 ~ ^ 7 ^ ~ ^ 2—6 ~ ^ ^ 

*> 

On the basis of previous experience, it is believed that the 
maximum safe operating temperature for air in contact with 
aluminum-alloy structures is between 300° and 400° F. Con- 
sequently, the air-temperature rise above ambient air from 
the heated-air source was assumed to be 300° F. If a 10° 
air-temperature drop is assumed between the heated-air source 
and station 82, the temperature of the air at the corrugation 
entrance will be 290° F. The air-temperature drop through 
the upper corrugation from 0- to 3-percent chord was found by 
trial and error to be 87° F. Therefore, the average tempera- 
ture of the air in the corrugation will be 

( t ) = 290 - — = 346° F 

2 av 2 

Since the design free-air temperature t 7 is 0° F, the re- 
quired skin temperature t 6 will he 100° F. The average 
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external heat-transfer coefficient h g _ ? f° r the 0** to 3- 

percent -chord upper surface, obtained by graphical integra- 
tion from figure 15, is 15.1 Btu per hour per square foot, - 
degree Fahrenheit. The required internal heat-transfer coef 
ficient, as obtained from equation (7), is as follows: 


. h 6 , 7 <t 6 - * 7 > 

2 ” 6 " t - t 

Z 2 Z 6 


( 8 ) 


Substituting the above values into equation (8) gives 


h 


3-6 


15.1 x (100 - 0) 
246 - 100 


10.34 Btu per hour, square 

foot, degree Fahrenheit 


The equivalent diameter for a corrugation is 0.0192 foot. 
Using the value of thermal conductivity for air at 246° F 
obtained from figure 10, the Nusselt number is 

Hu = h 3-s D ° = 10-34 . -* 0-0193 = 10 .30 
k 0.0193 


The corresponding Reynolds number is taken from ^figure 24, 
which gives the Reynolds and Nusselt number relationship ob- 
tained from figure 76, page 171, of reference 8 for the flow 
of air inside pipes. 

GD„ 3 

Re = — 2- = 2.50 x 10 3 

U- 

Using the value of absolute viscosity for air at 246° F as 
obtained from figure 10 gives 

3 151x1 0” 5 

G = 2.50 x 10 x — : = 1.965 pounds per second, 

0.0192 square foot 


The cross-sectional area of a single corrugation is 0.000651 
square foot. The weight-rate flov; of heated air required to 
raise the temperature of the strip of skin analyzed 100° F, 
then, is - 

w = G x A = 1.965 x 3600 x 0.000651 = 4.61 pounds per hour 

per corrugation 
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To check the heat balance, the total heat transferred from 
the strip of skin to the ambient air is assumed to be equal 
to the total heat removed from the heated air in passing 
through the corrugation from 0- to 3-percent chord: 

/ h 6 _ 7 S(t 6 - t 7 ) = wc p At a (9) 

where S = L/12 for a single corrugation. The value of L 
is obtained from figure 25. Then 

h c _ — (t R - t_) = 15.1 x 0 - x 100 =-97.4 Btu per hour 
e-? 12 6 7 12 


and 


wCp At 2 = 4.61 x 0.242 X 87 = 97.0 Btu per hour 


Similar calculations were made for the wing upper surface 
at station 382. The air-temperature drop from stations 82 to 
382 was assumed to be 30° F, which on the basis of flight data 
from similar designs appeared to be a reasonable value. Based 
on this assumption, the temperature of the air entering the 
corrugation then will be 260° F, The air flow necessary to 
raise the skin temperature 100° F for the 0— to 3-percent- 
chord section was found to be 7.01 pounds per hour per corru- 
gation. The analysis was completed for the upper corrugations 
of wing stations 82 and 382, using the previously established 
air-flow rates to determine the temperature rise of the skin 
aft of the 3-percent-chord point. 


The air-pressure drops through the upper corrugations at 
stations 82 and 382 were calculated for these air-flow rates, 
using the equation 


Pi " P 3 = 


+ fav Y av 


S 


2gm 


( 10 ) 


which is a form of equation (22c), page 128 of reference 8, 
where the subscripts l and 3 represent the regions shown 
in figures 3 and 5. The values of the friction factor f av 

were obtained from figure 26, which is taken from figure 51, 
page 118 of reference 8. It was found that the pressure loss 
through the corrugation at station 382 was greater than that 
at station 82, and therefore this air distribution never 
could be obtained with the same size of corrugations at the 
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tv/o stations. Consequently, the surface temperature for the 
section from 0- to 3-percent chord could not he raised ex- 
actly 100° F simultaneously at the two stations. Previous 
experience has shown that after an approximate required value 
of air-flow rate has been found, based on the design 'skin- 
temperature rise, the actual air-flow rate which will give 
pressure equilibrium is most readily established by assuming 
air-flow rates, checking pressure drops by an approximate 
method, and then calculating the skin-temperature rises. 

This procedure, therefore, was followed for the remaining 
wing- out er-panel calculations. 

Recalculation of heated-air distribution. - The approxi- 
mate method of determining the pressure loss in the corruga- 
tion passages consisted in using the second term only of equa- 
tion (10), since the first term is of minor importance and 
can be introduced as a refinement of the final calculations 
without affecting the results to any appreciable extent. 

For the upper corrugation at wing station 82 an air-flow 
rate, based on the preliminary calculations, of 5.00 pounds 
per hour was assumed. Then 

0 = — = — — — = 2.13 pounds uer second, square foot 

A 0.000651 X 3600 


Based on this air-flow rate, the air-temperature drop in the 
corrugation from 0- to 3-percent chord was found by trial to 
be 84° F. Then the average air temperature is 

(t 2 ) = 290 - — = 248° F 

3 av 2 

Using the value of absolute viscosity for air at 248° F, ob- 
tained from figure 10, the Reynolds number is 


Re 


£5e 


2.13 X 0.0192 
1.52 X 10^ 5 


2.69 X 10 3 


and from figure 24 


h__ D 

Nu = -2—^ = 10.9 

k 

Using the value of thermal conductivity at 248° F, obtained 
from figure 10, the available internal heat-transfer coeffi- 
cient is 
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■ h - = - — * - — - — — ' = 10.96 Btu per hour, square foot, 

° 0.0192 degree Fahrenheit 

From figure 15, 

h 6 _ 7 = 15.1 Btu per hour, square foot, degree Fahrenheit 


The average skin-temperature rise is calculated from the fol- 
lowing equation: 


* 6 “ 


h 3-6 *3 

^ 2-6 + ^ 6-7 


( 11 ) 


which is reduced from equation (?) when t = 0° F . Thus, 


10.96 X 248 
10.96 + 15.1 


104° F 


To check the heat balance , the total heat transferred from 
the 1-inch strip of skin to the ambient air is assumed to he 
eoual to the heat transferred from the heated air in travel- 
ing from 0- to 3-percent chord (equation (9)). Substituting 
appropriate values gives 

h _ _ — (t„ - t „ ) = 15.1 X 0 ' - -- 4 - X 104 = 101.4 Btu per hour 
e-? 12 6 7 12 


and 


wCp At 2 = 5.00 X 0.242 X 84 = 101.6 Btu per hour 

which show that thermal equilibrium has been established for 
this section. 

Repeating the calculations for the sections aft of 3- 
percent chord, using the values of external heat-transfer 
coefficient from figure 15, the skin-temperature rises shown 
in table III for the upper surface at station 82 were ob- 
tained . 

For an upper corrugation at station 382 the entering air 
temperature was assumed to be 260° F and the air-flow rate to 
be 6.20 pounds per hour per corrugation. The resulting skin- 
temperature rises were calculated and are presented in table 
III. 
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Pressure-loss calculations were made for the upper cor- 
rugations at stations 82 and 382 to determine if the pressure 
drop at station 382 was sufficiently less than the pressure, 
drop at station 82 to allow a reasonable pressure drop 
through the leading-edge duct. In the succeeding pressure- 
drop calculations it will be shown that the required pressure 
equilibrium is attained with these established air-flov; rates. 

Before final pressure-loss calculations could be made, 
it was necessary to determine the air-flow rates through the 
lower corrugations at stations 82 and 382 and through the 
corrugations at the intermediate stations. The lower-corru- 
gation air-flov; rates are dependent upon the upper-corrugation 
flow rates, since the air-pressure drops through the upper 
and lower corrugations must be equal. After determining the 
necessary 1 ower-c orrugat i on air-flow rates by the process of 
trial and error, the resulting skin-temperature rises were 
calculated by the method previously outlined, using the val- 
ues of L obtained from figure 27. 

Arbitrary curves of heated-air-flow distribution, the 
general shape of which was based on experience, were then 
drawn through the points representing the flow rates estab- 
lished for stations 82 and 382. The assumed curves of the 
flow distribution through the wing leading— edge corrugations 
are shown in figure 28. The correct shape of these curves is 
governed by the interdependency of the flow of heated air 
through the leading-edge duct and corrugations and the air- 
pressure drops through the leading-edge duct and corrugations, 
and will be discussed later. On the basis of these curves, 
skin-temperature rises were calculated for stations 157 and 
292, using the values of heat-transfer coefficients from fig- 
ures 16 and 17 . The calculated wing skin— temperature rises 
for all four wing stations are given in table" III. By inte- 
grating the areas under the heated-air-flow distribution 
curves, the total required flow of heated air was obtained. 

In determining the total heated-air-flow rate the fact that 
nose ribs would be spaced 15 inches apart in the wing, thereby 
eliminating a corrugation every 15 inches, was taken into 
consideration. There are appr oximat- ely 370 corrugations in 
the wing outer panel from station 11 to station 412. 

Calculate from figure 28 the weight rate of air flow 


^ a ^total = pounds per hour 
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Using the assumed air-temperature rise through the heat ex- 
changer (300 F), the required exchanger thermal output can 
he established as follows: 


^exch = ^a c p = ^130 x 0.243 X 300 = 301,000 Btu per hour 


Curves of the total heat removed from the wing leading- 
edge surface, obtained by plotting summations of the heat re- 
moved from each corrugation analyzed, are shown in figure 28. 
The total heat removed, obtained by integrating the areas 
under the curves, was found to be 115,500 Btu per hour, or 
38 percent of the total heat supplied. 

Air-pressure drops in corrugation passages .- The follow- 
ing air-pressure-drop calculations are the results of a trial 
and-error development which was made in conjunction with the 
previous air-flow and heat-flow calculations. 

The air-pressure drops in the corrugations were calcu- 
lated, using equation (10). The air-pressure drop through 
the upper corrugation at wing station 82 was found as follows 

The air-flow rate for the upper corrugation at station 
82, as determined previously, is 5.00 p ound s per hour and 
® = 2.13 pounds per second, square foot. The air temperature 
at the corrugation inlet was assumed to be 290° F. The air 
temperature at the corrugation outlet was computed in the 
heat— flow calculations for the 0- to 10-percent-chord section 
and was found to be 147° F. Thus 

= t x + 460 = 290 + 460 = 750° F absolute 

and 

T 3 = t 3 + 460 = 147 + 460 = 607° F absolute 


and 


Then 


av 


RTj _ 53.3 X 7 50 
P 0 1058~ 

RT 3 _ 53.3 x 607 
~V~ 0 1058 

Vj + v 3 ^ 37.8 + 
2 2 


37.8 cubic feet per pound 


30.6 cubic feet per pound 


30.6 

= 34.2 cubic feet per 

pound 
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The value of (—^) 

V U / av 

82 was assumed equal to 
region, 3 percent to 7 


for the upper 

the value of 
percent , which 


corrugat i on 
from table 


at station 

the central 
III is 


( ) = 2.88 X 10 3 

V M- /av 

From figure 26, 

f = 0.0114 

The value of N in the case of the corrugations is equal to 
L, and is taken from figure 25. For a single corrugation 
m = 0.0048 foot. Then the pressure drop from region 1 to 
region 3, neglecting entrance and exit losses, from equation 
(10) is 

p - p = (30.6 - 37.8) + 0,0114 X 1.893 x 2.13 3 X 34.2 . 

1 * 3 32.2 2 X 32.2 X 0.0048 

= -1.015 + 10.85 


= 9.84 pounds per square foot or 1.89 inches of water 


Similar calculations were made for the lower corrugation 
at station 82 with a flow rate of 5.20 pounds per hour per 
corrugation. The air-pressure drop through the corrugation 
was found to be 1.90 inches of water, which is substantially 
the same as the value for the upper corrugation, indicating 
that pressure equilibrium for the upper and lower corrugations 
has been established. 

Repeating the calculations for wing station 382, the . 
pressure drop through the upper corrugation was found to be 
1.32 inches of water with a flow rate of 6.20 pounds per hour, 
and 1*35 inches of water through the lower corrugation with a 
flow rate of 6.70 pounds per hour. A tabulation of the pres- 
sure-drop calculations for the four wing stations analyzed is 
presented in table IV. 

Air-pressure drop through wing leading-edge duct . - The 
air-pressure-drop calculations for the wing and empennage 
leading edges were based on the assumption that the air trav- 
eled^through a smooth duct formed by the baffle plate and the 
leading— edge corrugations. In the case of the wing, the 
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leading-edge duct was considered to be formed by the baffle 
plate and the nose-rib liner for the section of span inboard 
of station 292. 

Using figure 28, the average air-flow rate from station 
11 to station 82 for the wing is 5.00 pounds per hour per cor- 
rugation for the upper surface and 5.15 pounds per hour per 
corrugation for the lower surface, or a total average of 

pounds per hour per corrugation. Since there are five 
nose ribs between stations 11 and 82, the total air lost 
through the corrugations from station 11 to station 82 is 

10.15 x [(82 - ll) - 5] = 670 pounds per hour 

Then the air-flow rate in the leading-edge duct at station 82 
i s 

W a = 4130 - 670 = 3460 pounds per hour 


The vring 1 eading- edge-duct dimensions are given in figure 29. 
At station 82, 


A = 0.458 square foot 
m = 0.160 square, foot 


Then 


G- 


W a 3460 

A 3600 X 0.458 


Ug = 4ra = 0.64 foot 


2.10 pounds per second, 
square foot 


When 


tj. = 290° F, m- = 1.59 x 10 5 pound per second, foot 


and 


_ GD e 2.10 x 0.64 4 

Ee - IT ■ 1.58 X 10- = ’ 8 - 50 X 10 


From figure 26, 
f = 0.00483 


<*i> 


RT X _ 53 .3 X 750 


av p 0 


1058 


= 37.8 cubic feet per pound 
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To calculate the air-pressure drop, equation (lO) was used. 
Since the air-temperature change is small through the incre- 
ments of duct analyzed, the first term may he neglected. 

Then the pressure drop at station 82 for 1 foot of span is 

fNG 3 v , r 0.00483 X 1 X 2.10 s X 37.8 
AP X = Si = 

2gm 2 X 32.2 X 0.160 

= 0.0781 pound per square foot per foot or 0.01501 
inch of water per foot 


Repeating the calculations for several other stations, as 
shown in table IV, plotting the results (fig. 30), and inte- 
grating the area under the curve from stations 82 to 382 re- 
sults in a total air-pressure drop through the leading-edge 
duct between these stations of 0.52 inch of water. For equi- 
librium conditions, the air-pressure drop through the corru- 
gations at station 82 must be equal to the air-pressure drop 
through the leading-edge duct from stations 82 to 382 plus the 
air-pressure drop through the corrugations at station 382. 




3 82 


From the calculations, 

(p - p ) = 1.90 inches of water 

1 r 3 82 


( Ap T ) P3 t0 382 + (pi - 23)302 = °- 52 + 1,34 = 1-86 inches of 

water 

which show that pressure-equilibrium conditions have been 
adequately established. Since (p a - p 3 ) sa slightly 

greater than (Ap 1 ) ea tQ 383 + (p a - P 3 ) 382 . a slightly 

greater air flow through the corrugations near station 382 
might be expected, with a corresponding slight increase in 
skin-temperature rise in the region near the wing tip which, 
as seen from table III, is desirable. 

A comparison of the pressure-drop difference of the air 
passing through the corrugations at the four wing stations 
with the air-pressure drop through the leading-edge duct, with 
station 82 as reference, is shown in figure 31. The shape of 
the air-flow-distribution curves (fig.. 28) depends on these 
two curves. When pressure-drop equilibrium has been estab- 
lished, the two curves of figure 31 should coincide, since 
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the air-pressure-drop difference of any two corrugations must 
equal the air-oressure drop through the leading-edge duct be- 
tween those two points. From the curves shown in figure 31, it 
is evident that the air- f 1 ow-d i str i buti on curves should be 
slightly less concave in order to establish equilibrium. 

The method of calculating the air and heat flows and 
pressure drops for the horizontal stabilizer and vertical fin 
was the same as for the wing. A 1 e ad ing- edge- duct air tem- 
perature of 270° F was assumed at horizontal-stabilizer sta- 
tion 72 and at vertical-fin station 122. The general dimen- 
sions of the h or i z ont al- st ab i 1 i zer leading-edge thermal sys- 
tem as used in the calculations are given in figures 32 and 
33. The stabilizer air and heat flows and skin-temperature 
rises are presented in table V and figure 34, and the air- 
pressure drops are presented in table VI and figure 35. The 
general dimensions of the vertical fin as used in the calcu- 
lations are given in figures 36 and 37. The fin air and heat 
flows and skin-temperature rises are presented in table VII 
and figure 38, and the air-pressure drops are presented in 
table Vll and fi gu re 3 5. 


AIR-HEATED- WINDSHIELD DESIGN 


The approach to the air-heat ed-winfl shi e Id design was dif- 
ferent from that used in the analysis of the airfoil sections. 
Since the value of the external heat-transfer coefficient for . 
the wind shield surface was not known, the method outlined in 
reference 9 was used and a specified heat flow through the 
windshield outer panel was established as the heat requirement 
for ice prevention. 

As in the analysis of the wings and empennage, certain 
factors were fixed and an attempt was made to vary the remain- 
ing factors to satisfy the design requirements. The fixed 
factors consisted of the v/indshield area and shape, while the 
variable factors consisted of the heated-air temperature., the 
flow rate, and the gap thickness. Consideration of the 
strength of the glass and the comfort of the pilots limited 
the heated-air temperature to a maximum value of about 200° F . 

The necessary heat flow through the windshield outer 
panel was taken as 1000 Btu per hour per square foot, as rec- 
ommended in reference 9. The cuter panel consists of two 
layers of semi- 1 emper ed plate glass 0.1095 inch thick each, 
separated by a layer of vinyl plastic 0.125 inch thick. The 
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thermal conductivity of the glass was taken as 6.67 Btu per 
hour, square foot, degree Fahrenheit per inch and the thermal 
conductivity of the plastic was assumed to he 1.25 Btu per 
hour, square foot, degree Fahrenheit per inch, although this 
value is not considered accurately established. 'The thermal 
conductance C for the outer panel was calculated from equa- 
tion (3) of reference 9; 


0 = - ( 12 ) 

I i I 2 ^3 

r, + y; + k; 


where l designates the thickness and K the thermal con- 
ductivity of each of the three layers. Then 


C 


1 

0.1095 x 2 + 0.125 
6.67 1~ 25 


7.54 Btu per hour, square foot, 
degree Fahrenheit 


In order to obtain the required heat flow through the outer 
panel at the outer-surface temperature of 50° F, as specified 
in reference 9, the inner-surface temperat’.xre must be 


t 


i 


1000 4 50 0 

c 


(13) 


which for the subject windshield is 183^ F. As was stated 
previously, the heated-windshield design included a secondary 
heat exchanger which would transmit the heat of primary heat- 
ing air from the exhaust exchangers to secondary air for 
windshield and cockpit heating. Calculations using the method 
of reference 9 indicated that no practical double-panel design 
which would give the required heat flow was possible with the 
heat available from a secondary heat-exchanger installation 
that used the remaining heated air from the inboard exchangers 
after supplying the empennage system. The secondary exchanger 
installation was necessary, however, in order to supply cock- 
pit heated air free from carbon-monoxide contamination in the 
event of leakages in the exhaust exchangers. The windshield 
installation finally adopted consisted of (l) the double-panel 
arrangement which would provide the optimum heating for the 
allowable pressure drop through the gap, and (2) a means of 
discharging the primary air from the secondary exchanger over 
the outer surface of the windshield. In order to provide for 
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circulation of air through the secondary exchanger and wind- 
shield system by the dynamic- pressure of the free stream, 
the allowable pressure drop through the windshield gap was 4 
inches of water . Based on the equations of reference 9, it 
was indicated that the optimum heating of the windshield 
outer panel with a 4-inch pressure drop in the gap would be 
provided with an air-gap thickness of approximately 3/16 inch. 

The reason for the difficulty in obtaining the necessary 
heat flow through the outer panel was due in part to the lack 
of heat available and in part to the low thermal conductance 
of the panel. This low conductance was caused by the rela- 
tively thick plastic layer of low thermal conductivity. 
Therefore, it is apparent that in the design of air-heated 
windshields attention should be directed to securing an outer 
panel of high thermal conductivity. 


RESULTS AND ERRORS 


A general summary of the calculated heated-air flows, 
skin-temperature rises, and air-pressure drops for the wing, 
horizontal stabilizer, and vertical-fin thermal systems is 
presented in figures 40, 41, and 42. Table VIII gives a sum- 
mary of the calculated heat flows. 

Aside from experimental errors involved in the determi- 
nation of airfoil pressure distribution, the foregoing treat- 
ment of the analytical problem has introduced several errors 
which should be acknowledged. 

An error was introduced in the use of the arithmetic 
average for obtaining ^2^ a v* the avera £ e air temperature 

in a corrugation between two points. A calculation of the^ 
air-temperature change through the upper corrugation at wing 
station 82 was made using small increments of corrugation 
length. This more exact calculation of t 2 is compared in 

figure 43 with the straight-line temperature gradient assumed 
in the analysis. The maximum deviation between the two 
curves was found to be 1,24 percent. A tabulation of the 
calculations made to obtain the more exact curve of the air- 
temperature change is presented in table IX. 

The effects of radiation and conduction from the corru- 
gation walls to the leading-edge skin were also neglected. 
However, since the corrugation ducts are 3/4 inch wide and 
are required to heat a 1-inch strip of skin, these errors 
tend to compensate each other. 



NACA ARK No. 5A03 


28 


In the air-pressure-drop calculations, the corrugation 
entrance and exit losses were neglected. In the case of the 
wing, the pressure drop of the air traveling from the nose- 
rib-liner duct to the corrugation entrances was neglected as 
we 11. Since the purpose of the pressure-loss calculations 
was to determine the air-flow distribution which would result 
during flight, neglecting the entrance and exit losses would 
change the air- and heat-flow calculations only slightly. 

The air velocities in the corrugations at the tip regions are 
higher than in the corrugations at the root regions; conse- 
quently the entrance and exit losses would he higher in the 
tip regions than at the root sections. The air-flow rates 
through these corrugations, therefore, would probably be less 
than calculated. The curve of friction factor f as used in 
the calculations of pressure drop is intended for isothermal 
flow, as stated in reference 8, but was assumed to be suita- 
ble for the polytropic flows of the analysis. 


BLOWER TESTS 0? WING LEADING-EDGE THERMAL SYSTEM 


Upon completion of the right-wing leading-edge thermal 
system, blower tests of the leading edge were conducted to 
obtain isothermal air-flow data to determine the desirability 
of installing a nrrse^rib liner , 

A blower was connected to the leading-^dge duct with a 
venturi meter installed to obtain the air-flow rate. With 
appr oxi mat ely the required air-flow rate of the analysis sup- 
plied to the leading-edge system by the blower, the total 
pressure of the air out of representative upper and lower 
corrugations between each nose^-rib bay was measured by means 
of a survey tube connected to an open-end manometer. Prom 
these total pressures, the air-flow rates in the corrugations 
were calculated. 

A blower test was first made without t he nose-^r ib liner 
installed in order to determine if a liner would be necessary. 
It was found that the pressure drop along the leading-edge 
duct was excessive and prevented the desired flow of air 
through the corrugations in the region of station 382. The 
nose-^r ib liner was then installed and holes were cut in the 
liner to give appr oxi mat ely the same area as the corrugations 
they supplied. Tests were made and the ai r- f 1 ow-di st r i but i on 
curve then compared favorably with the calculated distribu- 
tion curve. A comparison of the curves of the two tests with 
the calculated curves is shown in figure 44. 
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Although the experimental distribution curve is for iso- 
thermal air flow and the calculated distribution is for a 
polytropic flow, the tests indicate that the desired air-fl.ow 
distribution will be approximated at the design flight condi- 
tions. 


DISCUSSION 


The analytic procedure involved the employment of 
assumed atmospheric and flight conditions previously estab- 
lished as criteria in former design calculations of thermal 
ice-prevention systems. Tests of the equipment for the C-46 
airplane will help to determine the completeness of the de- 
sign criteria and the accuracy of the analytic methods, and 
possibly will permit their extension and refinement. The 
calculation of heat transfer from the airfoil surfaces to dry 
air was an expedient necessitated by a paucity of data on the 
water content and drop sizes in the atmosphere at various 
temperatures. Efforts are being made to determine from sta- 
tistical researches of the atmosphere values of temperature, 
pressure, water content, drop size, and humidity, which may 
be used to establish the optimum design criteria for thermal 
ice prevention. When such data have been established, the 
heat required can be based on the combined transfer due to 
air convection and water-drop impingement on the airfoil sur— ' 
faces. 

The use of the method of utilizing the airfoil pressure 
distribution to obtain the required heat for thermal ice 
prevention appears entirely feasible in the design of future 
thermal ice-prevention systems. The pressure distribution 
for any airfoil may be calculated as well as determined exper- 
imentally. Methods of calculating the external heat-transfer 
coefficients in a turbulent boundary layer are available 
(references 10 and ll) in the event that the transition point 
from laminar to turbulent flow occurs within the heated region 
of the airfoil. Once the required heat flows based on a given 
surface-t emperature rise are established, rough calculations 
can be made for various thermal-system designs to determine 
the optimum required heated-air flows. After the design has 
been decided upon, a complete analysis can be made. 

Summarizing, the general procedure for analysis is as 
follows: 

1. Determine the external heat-transfer coefficients 
from airfoil pressure distributions. 
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2. Determine the required heat flows through the air- 

foil leading edges for a given surface-temperature 
rise. 

3. Determine the necessary heated-air flows and required 

heat-exchanger output, and correlate the air flows 
with the resulting air-pressure drops for a given 
thermal- syst em design. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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TABLE I.- EXPERIMENTAL AIRFOIL PRES8URE DATA OBTAINED BY MEANS 
OF PRESSURE BELTS FROM FLIGHT TESTS OF THE C-46 AIRPLANE 
[l-P « l-Pf/q©] 


Station 82 



Station 

157 

Station 292 



Station 

382 

q 0 - 11 

.3 In. H2O 

Qo = U- 

6 in. H2O 

Qo = 12. 

in. HgO 

Qo 

= 12.4 

in. H2O 

pressure altitude, 

preesure 

altitude, 

pressure 

altitude, 

pressure 

altitude , 

5000 ft 



5000 ft 

18,000 ft 



18 , 000 

ft 

Wing 

Fercent 



Percent 



Percent 




Percent 


1 V.. 

chord 


Upper 

cnord 


"*Upper 

chord 


'■“''Upper 

chord 


1 *upper 

11.8 


2.27 


14.1 


2.33 

14.2 


2.07 

19.3 


1.70 

7.4 


2.31 


9.1 


2.41 

9.5 


2.08 

13.1 


1.73 

4.7 


2.24 



5.9 


2.36 

6.3 


2.00 


8.8 


1.66 

3.3 


2.08 



4.3 


2.25 

4.3 


1.85 


6.0 


1.57 

1.6 


1.64 



2.3 


1.88 

2.7 


1.67 


3.9 


1.43 

1.0 


1.46 



1.5 


1.62 

1.6 


1.47 


2.4 


1.24 

.4 


1.02 


.7 


1.29 

.9 


1.13 


1.3 


1.00 

.1 


.58 



.3 


.90 

.4 


.79 


.6 


.75 



1 " p Low 




1 

“ p Lower 

2 




3 


.49 



er 











.03 


.08 



.04 


.08 

0 


.14 


0 


.06 

.6 


.27 



.5 


.29 

— 


i'^Lower 

- 

— 


l~ p Lower 

1.1 


.47 



1.3 


.54 

.2 


.26 

< 

.7 


.90 

2.1 


.68 


2.4 


.83 

.6 


.97 


1.4 


1.23 

3.3 


.35 


3.8 


.93 

1.4 


.89 


2.6 


■1.13 

4.75 


.98 



5.6 


1.03 

2.5 


1.22 


4.0 


1.19 

6.6 


1.03 


7.7 


1.10 

3.9 


1.16 


6.2 


1.21 

9.1 


1.06 


10.6 


1.15 

6.0 


1.25 


8.6 


1.11 

13.0 


1.16 


15.0 


1.21 

8.4 


1.22 

11.8 


1.15 









12.0 

17.0 


1.24 
1 so 

16.5 
s.v 0 


1.04 

1.01 


















** 





Station 72 

Station 127 


Station 172 


Station 122 

Station 172 

q 0 *= 9.6 in. H2O 

Qo = 

12. 

2 in. H2O 

Qo * 12. 

4 in. H2O 

Qo * 

3 in. 

H 3 0 

Qo * U. 

6 in. H2O 

pressure 

altitude, 

preesure 

altitude. 

pressure 

altitude, 

pressure 

altitude, 

pressure 

altitude, 

5000 

ft 


18,000 ft 


18,000 ft 


5000 ft 


5000 ft 

Horizontal stabilizer 

Vertical fin 

Percent 

chord 

!-Pupper 

Percent 

chord 

1_P Upper 

Percent 

chord 

l _p Upper 

Percent 

chord 

1 - p Left 

Percent 

chord 

1 - p Left 

13.8 


1.32 

18. 

4 

1.39 

22.5 

1.32 


7.11 

1 

29 

11.10 


1.35 

10.0 


1.30 

13. 

7 

1.42 

17.0 

1.10 


4.78 

1 

29 

7.49 


1.44 

6.3 


1.29 

9. 

1 

1.39 

11.4 

1.33 


3.34 

1 

24 

5.29 


1.43 

4.0 


1.23 

6. 

3 

1.37 

8.1 

1.29 


2.38 

1 

18 

3.82 


1.36 

2.5 


1.08 

4. 

4 

1.32 

5.8 

1.23 


1.46 

1 

04 

2.41 


1.27 

1.1 


.92 

2. 

6 

1.19 

3.6 

1.12 


1.03 


96 

1.69 


1.15 

.4 


.58 

1. 

7 

1.03 

2.4 

.99 


.81 


90 

1.36 


1.03 

.2 


.25 

1. 

3 

.93 

1.9 

.93 


.61 


82 

1.04 


.93 

.01 


.02 


8 

.76 

1.3 

.82 


.41 


72 

.64 


.74 


1 

“ p Lower 


4 

.44 

.7 

.61 


.11 


40 

.14 


.38 

.01 


.14 

A 


A 

U 


.03 

-.06 


0 


31 



i . . . 











w Rigbt 

.5 


.90 


- 

1-P Lower 


l -p Lower 


— 

1 - p Right 

.02 


.58 

1.2 


1.13 


5 

1.03 

.4 

1.06 


.03 


40 

.17 


.78 

1.6 


1.14 

1. 

0 

1.17 

1.0 

1.22 


.11 


53 

.42 


.88 

2.0 


1.17 

1. 

4 

1.32 

1.5 

1.36 


.23 


69 

.69 


.97 

2.4 


1.18 

1. 

9 

1.33 

2.1 

1.39 


.61 


94 

1.36 


1.15 

3.2 


1.25 

2. 

B 

1.40 , 

3.2 

1.43 


1.46 

1 

10 

2.74 


1.18 

4.7 


1.30 

4. 

6 

1.51 

5.6 

1.50 


2.38 

1 

20 

4.16 


1.25 

6.3 


1.34 

6. 

6 

1.60 

7.8 

1.50 


3.79 

1 

27 

6.30 


— 

7.9 


1.34 

9. 

3 

1.57 

11.1 

1.45 


6.16 

1 

33 

9.89 


1.27 

12.4 


1.31 

12. 

7 

1.51 

16.7 

1.43 


— 

— 

— 

— 


— 

16.2 


1.30 

17. 

4 

.1.55 

22.2 

1.41 


1 

— 




— 
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TABLE II.- DETERMINATION OF EXTERNAL HEAT-TRANSFER COEFFICIENTS 
FROM AIRFOIL PRESSURE-DISTRIBUTION MEASUREMENTS AT 
WING STATION 82 OF THE C-46 AIRPLANE 
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o 

r~* 

X 

£ 

“ 

o 

o 

rA 

X 

o 

K 

P. 

1 

rA 

o 

> 

> 

« 

O 

•4 

CD 

c 

Cj 

c 

c 

< 

< 

C 

3 

-1 

> 

A 

0 

D 

D 

3 

3 

3 

o 

CD 

TJ 

o 

O 

> 

> 

o 

r» 

00 O 
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r* 

GO 

r- 

e • 

o 

> 

> 

9 

^ p 
uD 

/< 

f — i 

c*. 

o 

p 

<M 

r- cr 
I CO 

c — 

S3 

u 

-d 

"5“ 

•p 

• i 

r — i 

u. 

o 

U P 
M <M 

0 cr 

. 03 

^ ' 

1 '~- 
C M 

S3 S3 

3 

p 

m 

i — i 

Upper surface 

18.88 

15.0 

2.196 

1.48 

8. 

550 

0.839 

2.675 

0.412 

2.00 

2.03 

13.77 

10.0 

2.303 

1.52 

6. 

070 

.598 

1.903 

.674 

4.60 

4.65 

11.69 

8.0 

2.312 

1.52 

5. 

155 

.511 

1.622 

.765 

6.13 

6.20 

9.54 

6.0 

2.291 

1 . 51 

4. 

240 

.423 

1.338 

.765 

7.44 

7.54 

7.29 

4.0 

2.185 

1.48 

3. 

310 

.324 

1.033 . 

.765 

9.62 

9.75 

6.11 

3.0 

2.030 

1.43 

2. 

880 

.317 

.955 

.765 

10.41 

10.55 

5.47 

2.5 

1.930 

1.39 

2. 

640 

.293 

.880 

. 765 

11.30 

11.44 

4.83 

2.0 

1.802 

1.34 

2 . 

415 

.265 

.829 

.765 

12.00 

12.15 

4.12 

1.5 

1.650 

1.29 

2. 

155 

.260 

.743 

.765 

13.30 

13.47 

3.37 

1.0 

1.450 

1.20 

1. 

875 

.241 

.672 

.765 

14.80 

15.00 

2.51 

.5 

1.140 

1.07 

1. 

570 

.233 

.605 

.765 

16.43 

16.64 

.80 

0 

. 250 

. 50 

— 

— 

— 

-- 


.765 



0 

.1 

0 

— 

-- 

— 

— 

-- 

.276 

.765 

36.00 

36.50 

Lower surface ”| 

0.99 

0.5 

0.230 

0.48 

1. 

390 


- 


0.765 



1.75 

1.0 

.425 

.65 

1. 

S06 

.375 

.822 

.765 

12.10 

12.27 

2.41 

1.5 

.560 

.75 

2. 

165 

.311 

.819 

.765 

12.14 

12.31 

3.02 

2.0 

.661 

.81 

2. 

500 

.251 

.792 

.765 

12.55 

12.71 

3.58 

2.5 

.743 

.86 

2. 

792 

.273 

.874 

.765 

11.39 

11.53 

4.14 

3.0 

.813 

.90 

3. 

084 

.280 

.925 

.765 

10.75 

10.90 

5.21 

4.0 

.921 

.96 

3. 

340 

.295 

1.036 

.765 

9.60 

9.73 

7.31 

6.0 

1.024 

1.01 

4. 

850 

.391 

1.378 

.765 

7.22 

7.31 

9.36 

8.0 

1.065 

1.03 

6.070 

. 455 

1.662 

.765 

5.97 

6.05 

11.39 

10.0 

1.077 

1.04 

7. 

36 

.538 

1.990 

.765 

5.00 

5.06 

16.45 

15.0 

1.160 

1.08 

10. 

27 

.533 

2.335 

.765 

4.26 

4.31 


TABLE III.- HEAT -FLOW AND SKIN-TEUPERATURE-RISE ANALYSIS FOR THE WING 
OUTER PANEL LEaDING-EDGE THERMAL SYSTEM OF THE C-46 AIRPLANE 


St at Ton 

Surface and 
percent chord 

'a 

o 

•H 
P 
flJ 
bo 
3 
u 
u 
o 
* o 

M 

0) 

p- 

M 

XJ 

.Q 

rH 

r - 1 
P 

<M 

o* 

GO 

o 

0) 

00 

S3 

rH 

l—l 

— El. 

a o 

p 

i — * 
p 

«o — - 
O ■ — ' 
rA O 
0) 
X CD 

rA 
l — 1 

« 

i 

o 

iH 

X 

a>| 

sr 

4) 

o 

c 

1 

a 

S3 

M 

r— i 
■P 

<M 

u. 

o 

p 

cm 

m cr 
CD 

U 

S3 

P 

P 

CQ 

uT 

o 

p 

® or 

1 CD 
a ' 
S3 ^ 
H 
S3 

\ 

3 

P 

CQ 

« — i 

O 

P 

<M 

► cr 
» © 
O w 
S3 ~ 
u 
S3 

p 

p 

m 
«_ < 

Average skin tempera- 
ture rise 
(°F) 

82 

U 

o 

0-3 

5.00 

2.13 

248 

1.52 

2.69 

10 

.9 

0.0193 

10.96 

15.10 

104 


p< 

p< 

3-7 

5.00 

2.13 

186 

1.42 

2.88 

11 

.5 

.0179 

10.71 

8.55 

104 



7-10 

5.00 

2.13 

156 

1.38 

2.96 

11 

.8 

.0171 

10.50 

5.75 

101 


U 

« 

0-3 

5.20 

2.22 

250 

1.52 

2.80 

11 

.3 

0.0194 

11.41 

15.60 

106 


> 

o 

3-7 

5.20 

2.22 

192 

1.43 

2.98 

11 

.8 

.0180 

11.07 

8.53 

108 


tA 

7-10 

5.20 

2.22 

164 

1.39 

3.06 

12 

.1 

.0173 

10.90 

5.70 

108 

KS1 

U 

0-3 

5.05 

2.16 

242 

1.51 

2.75 

11 

.1 

0.0192 

11.10 

18.60 

91 

■ 

1 

3-7 

5.05 

2.16 

185 

1.42 

2.92 

11 

.6 

.0178 

10.74 

9.52 

98 

■ 

o 

7-10 

5.05 

2.16 

157 

1.38 

3.00 

11 

.9 

.0172 

10.67 

6.55 

97 

HI 

u 

0-3 

5.30 

2.26 

246 

1.51 

2.87 

11 

.5 

0.0193 

11.57 

16.60 

101 

. 

> 

3-7 

5.30 

2.26 

196 

1.44 

3.02 

12 

.0 

.0181 

11.31 

8.80 

110 


•H 

7-10 

5.30 

2.26 

170 

1.40 

3.10 

12 

.2 

.0175 

11.11 

6.42 

108 

292 

u 

0-3 

5.50 

2.35 

‘ 243 

1.51 

2.99 

11 

Tif 

0.0192 

11.90 

19.10 

93 


£ 

3-7 

5.50 

2.35 

202 

1.45 

3.11 

12 

.2 

.0182 

11.58 

10.43 

106 


& 

7-10 

5.50 

2.35 

179 

1.41 

3.20 

12 

.5 

.0177 

11.51 

7.78 

107 


u 

0-3 

5.85 

2.50 

246 

1.51 

3.18 

12 

.4 

0.0193 

12.48 

21.60 

90 


> 

3-7 

5.85 

2.50 

209 

1.46 

3.29 

12 

.7 

.0184 

12.17 

9.75 

116 


3 

7-10 

5.85 

2.50 

189 

1.43 

3.36 

13 

.0 

.0179 

12.11 

6.61 

122 

382 

U 

© 

0-3 

6.20 

2.65 

241 

1.51 

3.37 

13 

.0 

0.0192 

13.00 

20.50 

94 


a 

3-7 

6.20 

2.65 

211 

1.46 

3.48 

13.4 

.0184 

12.83 

11.84 

110 



7-10 

6.20 

2.65 

192 

1.43 

3.56 

13 

.6 

.0189 

12.74 

8.70 

114 


u 

© 

0-3 

6.70 

2.36 

244 

1.51 

3.64 

13 

.9 

0.0192 

13.90 

23.00 

92 


» 

3-7 

6.70 

2.86 

218 

1.47 

3.74 

14 

.2 

.0186 

13.75 

9.98 

126 


*3 

7-10 

6.70 

2.86 

203 

1.45 

3.78 

14 

.3 

.0183 

13.62 

5.10 

148 









TABLE IV.- HEATED-AIR PRESSURE-DROP ANALYSIS FOR WING OUTER PANEL CORRUGATIONS 
AND LEADING-EDGE DUCT OF THE C-46 AIRPLANE 
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TABLE V.- HEAT-FLOW AND SKIN-TEMPERATURE-RISE ANALYSIS FOR HORIZONTAL 
STABILIZER LEADING-EDGE THERMAL SYSTEM OF THE C-46 AIRPLANE 


NACA ARR No. 5A03 
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TABLE VI.- HEATED-AIR PRESSURE-DROP ANALYSIS FOR HORIZONTAL STABILIZER 
CORRUGATIONS AND LEADING-EDGE DUCT OF THE C-46 AIRPLANE 


NACA ARR No. 5A03 
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TABLE VIII.- SUMMARY OF HEAT-FLOW CALCULATIONS FOR THE 
WING, HORIZONTAL STABILIZER , AND VERTICAL FIN 
THERMAL SYSTEMS OF THE C-46 AIRPLANE 


Leading- 

edge 

surface 

w a 

(lb/hr) 

At 

air tem- 
perature 
rise 
(°F) 
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Total heat 
supplied to 
leading-edge 
thermal system 
above 0° F 
(Btu/hr) 

*R 

Total heat 
removed from 
leading-edge 
surface 
(Btu/hr) 

Wing 

Stabiliser 

Fin 

4130 

2514 

1120 

300 

275 

270 

301.000 

167.000 
73,000 

115,500 

32,640 

19,860 

Leading- 

edge 

surface 

Se x 100 
Q 

Percent beat 
removed from 
leading edge 

S 

Total heated 
leading-edge 
surface 
( sq ft) 

« R /S 

Average heat removed 
from leading-edge 
surface 

[Btu/(hr)(sq ft)] 

Wing 

Stabilizer 

Fin 

38 

20 

27 

94.5 

18.9 

14.7 

1220 

1730 

1350 


TABLE IX.- CALCULATION OF HEATED-AIR TEMPERATURE CHANGE 
THROUGH UPPER CORRUGATION AT WING STATION 82 
OF THE C-46 AIRPLANE 


w = 5.0 pounds per hour per corrugation 

G *= 2.13 pounds per second, square foot 

c = 14.51 feet 
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Fig. 1 



Figure 1.- The C-46 airplane In which the thermal 
Ice-prevention equipment was installed. 
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Fig. 7 
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Figure 8.- Pressure belt mounted on the C-46 vertical 

fin leading edge for measuring the pressure distribution 
over the leading edge of the fin. 
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Fig. 9 



Figure 9.- The water manometer used for Indicating pressures 
at the empennage pressure belts, with a typical pressure 
indication, C-56 airplane. 
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0 2 4 6 8 10 12 14 16 

x/ c x 100 

Figure 11.- The airfoil pressure distribution for wing station 82. C-46 airplane. 






Lower s/c x 100 Upper 

Figure 12.- Ratio of local air velocity to free-stream velocity and boundary-layer thickness for flow over wing 
station 82, C-46 airplane. 






Maximum 



Figure 13.- Ratio of the local velocity at laminar separation, Us» to the maximum local velocity 
Ui, as a function of the slopes of the double-roof velocity profile. (Values 
obtained from reference 6. ) C-46 airplane. 












data for wing station 82, C-46 airplane. 












Lower surface s/c x 100 Upper surface 

Figure 17.- Calculated external heat-transfer coefficient from pressure-belt 
data for wing station 2^2, C-46 airplane. 
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Lower surface s/c x 100 Upper surface 

Figure 18.- Calculated external heat-transfer coefficient from pressure-belt 
data for wing station 382, C-46 airplane. 



NACA ABB No. 5A03 


net 19 



data for horizontal stabilizer station 72, C-46 airplane* 





data for horizontal-stabilizer station 127, C-46 airplane. 




data' for horizontal- stabilizer station 172, C-46 airplane. 





Btu/hr, ft 



data for vertical -fin station 122, C-46 airplane. 
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8/c x 100 

Figure 23.- Calculated external heat -transfer coefficient from pres sure -fcelV 
data for vertical-fin station 172, C-46 airplane. 
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Figure 25.- Distance around leading edge measured chordwise from 0 percent chord point for wing upper 
face, C-46 airplane. 
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Figure 30.- Heated- air-pressure drop through wing leading-edge duct, C-46 airplane. 
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Wing station 

Figure 31.- Comparison of air-pressure drop through corrugations with air-pressure drop through 
wing leading-edge duct, C-46 airplane. 




Figure 32.- Distance around, leading edge measured chordwise from 0 percent chord point for horizontal- 
stabilizer upper and lower surfaces,, C-46 airplane. 
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Horisontal-Btabilizer station 
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Figure 37.- Vertical-fin leading-edge -duct area and hydraulic radius, C-46 airplane 
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-fin leading-edge duct, C-46 airplane. 
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Fig. 43 



Figure 43.- A comparison of the straight-line interpolation of the 

heated-air-temperature change as used in the analysis with 
a more exact calculation for air in the upper corrugation at wing sta- 
tion^, C-46 airplane. 
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Figure 44.- Comparison of analytical heated-air distribution through corrugations with isothermal 
blower tests of wing leading-edge thermal system with and without noserib liner, C-46 

airplane. 



